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Abstract. A large number of energetic electrons are generated during solar flares. They 
carry a substantial part of the flare released energy but how these electrons are created 
is not fully understood yet. 

This paper suggests that plasma motion in an active region in the photosphere is the 
source of large electric currents. These currents can be described by macroscopic circuits. 
Under special circumstances currents can establish in the corona along magnetic field 
lines. The energy released by these currents when moderate assumptions for the local 
conditions are made, is found be comparable to the flare energy. 
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1. Introduction 

In the course of solar flares a large amount of energy is suddenly released 
and transferred into local heating of the coronal plasma, mass motions (e.g., 
jets and coronal mass ejections), enhanced emission of both electromagnetic 
radiation (from the radio- up to the 7-ray range) and energetic particles 
(i.e., electrons, protons, and heavy ions). Energetic electrons are responsi- 
ble for the non-thermal radio and X-ray emission of the Sun. Moreover they 
carry a substantial part of the released flare energy (ILin and Hudson! . Il97ll : 
Lml Il974l : lEmslie et all 120041 ). Oberservations by RHESS$\ indicate that 
flare produced high energetic (> 20keV) electron fluxes F e are of the or- 



der of F P sa 10 36 V s related to a pow e r P P o f abo ut P P as 10 22 W = 10 29 ^g/ s 



((Smith and Smith! . Il96.4 Irlevvaertsl . Il974l : Ehl 1 19741 : lEmslie et all 12004 



1 RHESSI stands for the Ramaty High Energy Solar Spectroscopic Imager 
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Figure 1: The figure shows a flare ribbon as observed by TRACE spacecraft on Octo- 
ber 28, 2003. The contour plots are obtained by RHESSI, i.e., the loop-top soft X-ray 
source is located in the middle and the footpoint hard X-ray sources are located beside 
the loop-top source. 



Warmuth et all 120071 ). However it is still unclear how so many electrons 
are accelerated up to high energies within fractions of a second. 

In this paper a realistic flare scenario is theoretically modelled. It con- 
siders genuine parameters basing on average observations of flaring active 
regions, e.g., see the event on October 28, 2003 as presented by the TiMCiill 
image in Fig. [TJ This image also contains the RHESSI X-ray contour plots, 
i.e., the loop-top soft X-ray source is located in the middle of the picture, 
whereas the footpoint hard X-ray sources are located next to the loop-top 
source. The following conclusions can be drawn from this figure: The hard 
X-ray sources (Fig. [1]) have a diameter of about S = 10 x 10 6 m resulting in 
a source area A s = tt (<2>s/2) 2 = 78.5 x 10 12 m 2 . Both hard X-ray sources are 
separated from one another by about L s = 70 Mm (Fig. [p. If it is assumed 
that both hard X-ray footpoints are located at the same height above the 
photosphere and belong to one circular magnetic field line, then this mag- 
netic loop has an arc length of L co = ( nL s)/2 ~ 110 Mm. 

The average kinetic energy W of one energetic electron can be estimated 
by W « p e/F e = 1022 w /(i0 36 s _1 ) ~ 62.4 keV corresponding to an average ve- 



1 TRACE stands for the Transition Region and Coronal Explorer 
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Figure 2: Simplified sketch of the geometrical configuration right even at the flare ignition 
in the solar corona. The current density j is established in the corona as a result of charge 
separation in the photosphere. 



locity V = 0.454 c 136 Mm /s, where c represents the speed of light. This 
velocity can be used to retrieve the electron density iV acc of the acceler- 
ated electrons N acc F e /(2A s v) = 1.17 x 10 13 m~ 3 . The 2 in this last equa- 
tion originates from the fact that RHESSI usually observes two hard X-ray 
sources at the footpoints (as it can be seen in Fig.Q]). By assu ming a typica l 



electr on density N co = 10 m A in the flare region (see e.g.. lAschwandenl . 
20021 ). i.e., considering the density corresponding to the electron plasma 



frequency of about 300 MHz, it can be seen that only a fraction of the 
available electrons is finally accelerated, i.e., iV acc « 1.2%iV C0 . However the 
energy contained in the accelerated electrons in comparison to the plasma 
electrons' thermal energy (three degrees of freedom are assumed) in the 
flare region can be determined by (^accW)/(i.5Ar co fe B r) ~ 14.1%. The quan- 
tity &b stands for Boltzmann's constant. Here a typical flare temperature 
T = 40 MK has been adopted , which is a value obt ained from the photon 



fluxes observed by RHESSI (I Warmuth et all 120071 ). Of course these val 



ues represent only rough estimates, but they impressively demonstrate that 
during flares just a small number of electrons are accelerated to relatively 
high energies. 

Currently several different electron acceleration mechanisms in the solar 
corona are known. All of these mechanisms have the principle of acceler- 
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(a) Translation of Fig. [2] into a circuit (b) The simplified electric circuit is ex- 



diagram. 
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Figure 3: Electric circuit diagrams 



ation due to electric fields in common, but differ in the processes leading 
to the generation of the electric field. In the present paper the generation 
of a large scale DC electric field is discussed in terms of electric circuits, 



tion fe.e.. Alfven 


and Carlcivist. 1967: Sen and White. 1972: Martres et al. 


1973; 


Hevvaerts . 


1974: Obavashi. 


97,4 lAkasofu. 19791: Kan et all 


1983 


Melrose!. 19971: Zaitsev et al.. 


1998; 


Yang et al. 


. 2004: Zaitsev. 2005) 


Mo- 



tivated by these papers, the electric currents are investigated in order to 
obtain a mechanism for the acceleration of electrons to high energies. The 
basic idea of this mechanism is to generate the flare energy by photospheric 
plasma motion in active regions. This is in contradiction to the reconnection 
model in which the magnetic field energy in the corona is taken for the flare. 



2. Description of the Model 

Consider a bipolar active region in the photosphere. The magnetic loops 
built arc like structures connecting regions of different magnetic polarities. 
These regions are separated from each other by the magnetic neutral line. 
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Since the temperature in the photosphere is about 5.8 kK the plasma is 
only partly ionised, whereas in the corona the temperature is sufficiently 
high to ionise most of the elements completely. The active regions are re- 
lated to hig h magnetic fi eld c oncentration penetrating the photosphere. 
For instance Yang et al. ( 20041 ) have reported magnetic flux densities in 
the case of NOAA AR 10486 (October 29, 2003) reaching up to 0.15 T. 
They also mention photospheric flow motio n with velocities up to 1.6 km /s . 



Kitahara and Kurokawa . 1990l ; 



These observation s are n ot unusual (see e.g. _ 
Wang et all 1200.4 l2004h . Hence it is kn own that the photosph e ric m atter 
motion (e.g., the Evershed motion ( Evershedl . 1909 : Chitre . 1968T). cv- 



is m 



clonic (also called vor tex like) plasma motion (see e.g., Martres et al. . 19731 : 



Hofmann et all Il992h) and often high p l asma shear velocities can be seen i n 



the photosphere (jSchleicher et all 120031 : lYang et all 120041 : IXu et all hooi ) . 



In Fig. [2] a bipolar active region is schematically presented. The plane 
represents an area on the solar photosphere, where a bipolar active region 
is located. The black dotted arrows represent the magnetic flux density B 
connecting the regions of different magnetic polarities which are separated 
by the magnetic neutral line. Due to the photospheric plasma motion the 
Lorentz force qui Q d X B acts on the charges q of the plasma and leads to 
generation of an electric current. If the direction of the plasma velocity 
■Ujnd is reversed at the velocity separatrix (dashed line) in Fig. [2j as it is 
indicated by the velocity arrows u; n d) then the Lorentz force points either 
toward the velocity separatrix or away from it, depending on the directions 
of the magnetic flux vector and the plasma velocity vector. 

The generated current leads to electric circuits: Electric currents al- 
ways choose the path of the lowest resistance. Since the plasma resistivity 
is highly dependent on the temperature, the conductivity i n the corona is 
about 1090 times higher than in the photospheric plasma ( Spitzer . 19651 ). 
In such a plasma the charged particles propagate along the magnetic field 
lines corresponding to electric wires. If there is a magnetic connection be- 
tween two oppositely charged areas through the corona, possibly as a result 
of magnetic reconnect i on, an electric current can close the electric circuit 
( Alfven and Carlqvist . 1967 ; Heyvaerts . 19741 ) through the corona (Fig. [2]). 
Then an electric field occurs along the coronal magnetic field lines and acts 
on the electrons within the coronal loop, and accelerates them along the 
magnetic field up to high energies beyond 100 keV. 
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3. The Electric Circuit Model 

The model (see Fig. [2j) described in the previous section is translated into a 



system of electric circuits as drawn in Fig. 3(a) There are two electric power 
supplies u\ and u 2 for the both different regions of magnetic polarity at the 
bipolar active region. Each of them has its own internal resistor, namely 
r 3 and r 4 . The induced current can be closed via the photosphere of each 
region, i.e., via the resistors r\ and r 2 , and/or by an interconnection between 
these both regions, i.e., via the resistors r^i and r\ t 2- These interconnections 
can be established by both through the photosphere, i.e., via the resistors 
r§ and r%, and through the corona via the resistors rs and rj. The latter 
can only happen, if there is a magnetic connection present between both 
differently polarised regions of the active region through the corona. For 
simplicity, the resistors rs and r%, as well as rj and r§ are combined to 
r i,i = r5r6 /( r 5+ r 6) and r\ t 2 = r ^ r s/{r 7 +r s ). Note that the resistors r§ and 
are coronal resistors, whereas the other ones are located in the photosphere. 
Since the resistivity is much lower in the corona than in the photosphere, the 
relationship r n <C r a for all n £ {5, 7} and all o € {1, 2, 3, 4, 6, 8} is satisfied. 
Therefore r- h i and rj )2 become either r- h i « r§ or r- h i « tq and rj |2 ~ rj or 
H,2 ~ r 8 depending on whether there is a magnetic connection through the 
corona or not, respectively. 



Applying Kirchhoff's laws on the circuit presented in Fig. 3(a) the im- 
portant result is = iq is found, i.e., there are always two equal but oppositely 
directed currents connecting the circuits of both magnetic regions. An ex- 



tended evaluation of the circuit system (Fig. 3(a)) provides the value for 
the current 

^ _ n(r 2 + r 4 )ui - r 2 (ri + r 3 )u 2 

rm (r 2 + r 4 ) + r 2 r 4 (n + r 3 ) + (r^i + r; j2 ) (n + ^3) (^2 + r A ) 



interconnecting both regions of the active region (IOnell . l2f)08h . According to 



Eq. JT]) a fully symmetrical circuit, i.e., U\ = « 2 , r\ = r 2 , and r 3 = r 4 would 
directly lead to i§ = iq = 0. This means the electrical circuit would be com- 
pletely closed through the photosphere and no current would flow through 
the interconnecting resistors, neither through the coronal part (rs and r?), 
nor through the photospheric part (re and r$). However a minor asymmetry 
(e.g., caused by different plasma flow velocities in the photosphere u\ 7^ u 2 
and/or different values of the resistors r%, r 2 , r 3 , and r 4 ) would lead to the 
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(a) The absolute value of the Dreicer field 
| -Ed I is shown in dependence on the plasma 
temperature T and electron density N e for 

the Case /3d = Aherm- 



Figure 4'- Dreicer field and Spitzer 
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(b) The electric resistivity r\ is plotted in 
dependence on the temperature for two 
different electron densities. The dots in 
the diagram mark the conditions for the 
photosphere (N e = 4 x 10 19 m~ 3 ) and the 
corona (jV e = iV co = 10 15 m" 3 ). 

stivity both according to lOnell (2008) 



occurrence of such oppositely directed (Fig. 3(a) ) currents 15 = iq / inter- 
connecting both parts of the bipolar active region. The electrons building 
up these two currents can be understood as the source of the hard X-ray 
footpoint sources, which are usually observed by RHESSI (see Fig. [T]). The 
hard X-ray emis sion is considered to be generated by energetic electrons via 
bremsstrahlung ( Brown . 1971 . 19721 ). 

Since electric currents always choose the path of the lowest resistance, 
which is through the corona in the considered case, the complete electric 
circuit can be simplified to a single mesh as drawn in Fig. 3(b) In that 
circuit a DC power supply U is serially connected with four macroscopic 
resistors (R p h,i, Rph,2, and n^)- The two power supplies from Fig. |3(a 



are considered to be merged into the one power supply of Fig. |3(b)[ i.e 
U = u\ — U2- Since the coronal contributions of the interconnection resistors 



and r; 2 from Fig. 3(a) are much smaller than the photospheric ones, 
n,i « R co ,i and r- h2 « R co ,2 is obtained. For simplicity R ph = i? ph ,i = R p h,2 



and R co = R co 1 = R co 2 are chosen. 



A macroscopic resistor R = (v l )/a is given by its cross sectional area 
A, length L, and electric resistivity 77. The first two parameters (^4 and 
L) are questions of geometry, whereas the electric resistivity is a plasma 
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parameter depending strongly on the plasma temperature T and weakly 
on the electron density iV e (see Fig. |4(b)[ ). To choose these parameters, 
the example shown in Fig. [T] is employed: For instance the hard X-ray 
source's diameter S is assumed to be about S 10 x 10 6 m, whereas 
its depth d a is (according to the depth of the photosphere) considered 
to be d s 500km (see e.g., IPriestl . 1 19841 ). thus A ph = d s s = 5 x 10 12 m 2 



is obtained. L p h = 40 x 10 6 m is used for the length of the photospheric 
resistor and the distance between the two hard X-ray sources is about 
L s = 70 x 10 6 m (see Fig. Q]). Assuming the overlying magnetic loop to be 
a semicircle, the length L co = [nL$)/2 ~ 110 x 10 6 m is found. The cross sec- 
tional area of the loop can be estimated by ^4 S which according to Sec. [L] can 
be obtained by RHESSI observations, i.e., A s m 80 x 10 12 m 2 = A co . Ac- 
cording to Fig |4(b)| the electric resistivity in the photosphere and corona for 
these given parameters is r] p h = 9.13 x 10 -3 fim and i] co = 8.37 x 10~ 6 fim, 
respectively. With these parameters the values 

Rph = !M^L = 73.0 x 10- 9 Q (2) 
Aph 

Rco = H^l = 1L5xW -^ n ( 3 ) 



are obtained for the resistors of the circuit shown in Fig. 3(b) 
The induced voltage 

U = u iQd BL ph (4) 

is determined using Faraday's induction law, if u; nc j and B denote the speed 
of the photospheric flow and the local magnetic flux density, respectively. 
Kirchhoff's law provides 

/mesh ( H ) " 2^ + ^) ^2lV " Uind ~ 



ph Vph 



for the current (see Fig. 3(b)). Note that the value of the current is inde- 
pendent of the length of the photospheric resistor. If it is assumed that the 
flare power P e = 10 22 W (see Sec. [Q) is equal to the electric power in the 



coronal resistors of the circuit in Fig. 3(b)[ then 



W(II) = ^/^«2.08xl0 16 A (6) 
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can be derived by using U co = -R C o^mesh(ii)i where U co represents the po- 
tential drop in one of the coronal resistors. Since there are two resistors in 
the corona, the 2 appears in the denominator of the middle-part of Eq. ((H) 
and in the following denominator of the electric current estimation. Note 
that the current is in good agreement with the electric current of about 
(F e e)/2 ~ 8 x 10 16 A, which is generated by the observed energetic (> 20keV) 
electron flux of about F e « 10 36 V s (see Sec. [TJ). The quantity e stands for 
the elementary charge. By using Eqs. j2|), JJJ, J5)), and J6|) a constraint for 
the product U[ nd B can be found 

u ind B = ^fflL « 76.1 v/ m . (7) 

This requirement (Eq. (fTjl) can be fulfilled, e.g., for u; nc i~ 870 m /s and 
B Ri 0.087 T whi ch are reasonable conditions for the photosphere (see e.g., 
Yang et all 120041 ). Finally the electric field along the coronal loop is esti- 



mated to be E = — u ™/l co as -2.18 x 10~ 3 v / m . The potential drop at one 
of the coronal resistors is U co 240 kV, and Eq. (|4|) gives U ~ 3 x 10 9 V for 
the power supply's voltage. 

In order to asses whether the photospheric flow can support the flare the 
following estimation is done: The continuous photospheric plasma motion 
along the way dX = n; nc jdt in the time interval dt builds up an energy of 

dW s = {N ph , e A ph L ph )(eu- md B)dX, (8) 

due to the action of the Lorentz force eu- m A B. Here iVph, e , and iVph^ph-kph 
stand for the total electron number density in the photosphere, and the total 
number of electrons in the volume of x L p h, respectively. Therefore the 
power of the photospheric motion 

Pph = ^jj^ = (N phfi A ph L ph ) (eu ind B) u ind 

= eN ph;e A ph L ph B uf nd w 8.4 x 10 25 W (9) 

(when L ph = 4 x 10 7 m, u ind = dx/ dt = 870 m/ s , and B = 0.087 T, as intro- 
duced in Sec. [37] and A^ p h, e = 4 x 10 19 m -3 are used) is much higher than the 
required flare released power in the corona. Hence the photospheric motion 
possesses more than enough power needed for the flare. 
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Figure 5: A to unity normalised one dimensional relativistic Maxwellian velocity distri- 
bution / re i plotted for a temperature of T = 1.4 MK is presented as a function of the 
normalised electron velocity 0. The runaway region /3d < /3 is filled, i.e., about 2.3% 
of the total electrons are accelerated when an electric field of E » —2.18 x 10 -3 v /m is 
chosen. 



4. Electron Acceleration 

Electrons gain energy only by electric field acceleration. In this paper the 
electron acceleration is caused by the electric field related to the potential 
drop at the coronal resistor. Since the electric field is aligned along the 
magnetic field in the coronal loop, the equation of motion for an electrorJl 
can be written in its one-dimensional form 

dp i I 

— = — eE — m e sign[p] \D\. (10) 

The quantities t and p = {™>* c $)l^f\^P denote the time and the electron's 
momentum, respectively. The normalised electron velocity f3 is given by 
(3 = V/c with V as the electron's velocity. The qua ntity D stands fo r the 
electron's deceleration due to Coulomb collisions ( Onel et al. . 200?! ). By 



inserting the expression for the electron's momentum into Eq. (|10l) 

^jj = m e c 7 3 ^ = ~eE - m e sign[/5] \D\ (11) 



3 The electron is supposed to carry the negative elementary charge — e and to possess 
the rest mass m e . 
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is obtained, if 7 = i/\/i-/3 2 for the Lorentz factor 7 is used. The momen- 
tum change per arc length x along the magnetic field can be derived from 
Eq. dm), i.e., d v/dx = m e cy 3 ■ d P/dx = 1 /(5c ■ dp/dt. 

An electron travelling through an electron-proton plasma experiences 
Coulomb collisions, which always make the electron lose energy. Thus the 
electron momentum changes by deceleration, as described in the very last 
sum of Eq. (fTTj) . The Coulomb deceleration D in an electron-proton plasma 
has two contributions, namely the electron-electron interaction D e and the 
electron-proton interaction Z? p , i.e., D = D e + D p . Each of these contribu- 
tions is given by 



VZ € {e,p} : Di 



4vre2 (i/ me + i/ mi y 2 c 2 pf 



for A / (12) 



(see e.g.. lOnel et all 120071 ). The quantity eo stands for the permittivity of 
free space, whereas Z\ represents the charge numbeiB, mi the rest mass, 
and N[ the number density of the particle species I G {e for "electron", p for 
"proton"}. Pi represents the relative velocity of the electron with respect to 
the electrons and protons of the plasma, in which the electron propagates. 

Hence f3 e = \J ft 2 + 3/? t 2 herm and (3 P j3 are the relative velocities of a moving 
electron with respect to electrons and protons of the (background) plasma. 
Aherm = ^ ( k B T )/ (m e c 2 ) is the thermal velocity normalised to the speed of 

light. The Coulomb logarithm In [Aj] = In J ( a d+ 6 o,j)/(26§ ~) contains the 

Debye length Ad = \J^ k ^ T )/{N e e 2 ) and the Coulomb collision impact pa- 

rameter 6o,j 



^oc 2 (i /me+ i/ m .y 

As i t can be se en from Eq. 



Dreicerl (Il95fll . Il96(ll ) field 



Til a special electric field, the so called 



£ D =-(— signal 
V e 1 



(13) 



exists, which makes the time derivative of the electron's momentum vanish 
by definition, i.e., A v/dt = 0. Thus there is a critical velocity, the so-called 



In a fully ionised electron-proton plasma Z e — 1 and Z p = 1 is satisfied. 
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Dreicer velocity /3d, which characterises the change of the electron's mo- 
mentum: If the initial electron velocity |/3j is greater than the Dreicer ve- 
locity (|/3;| > /?d ), the electrons a re accele r ated nearly collisionless and 



are called runaway electrons (e.g., Holman . 19951 ). In the opposite case 



< |/ 3p|) they are thermalised due to the strong Coulomb collisions. 
The Fig. 4(a) shows the dependence of the absolute value of Ejy on N e and 
T for (3r> = Aherm- 

The previously determined electric field of E ~ —2.18 x 1CP 3 v /m in the 
corona corresponds to /3d ~ 2.83/3 t herm and therefore all the electrons which 
possess a normalised velocity higher than 2.83 Aherm are located in the run- 
away regime. If an one dimensional relativistic Maxwellian velocity distri- 
bution at a temperature of 1.4 MK is assumed then about 2.3% of the total 
available electrons are frictionless accelerated by this electric field. 



5. Conclusions 

Electron acceleration is a very important process during solar flares, since 
a large fr action of the released f l are energy is depos ited in the energetic 
electrons ( Lin and Hudson . 1971 : Emslie et al. . 20041 ). Such electrons can 



be observed by in-situ spacecraft measurements, by their non-thermal ra- 
dio signatures, and by the X-ray radiation which they emit, if they travel 
through the dense chromosphere. The paper at hand explains how a large 
scale DC electric field is generated in the corona, which subsequently can 
accelerate the electrons: The photospheric motion in bipolar active regions 
can generate electric powers high enough to drive a system of electric cir- 
cuits in the solar atmosphere. The photospheric motion can induce voltages, 
because the photospheric plasma is not fully ionised, i.e., it contains a sub- 
stantial part of neutral particles. Due to the low resistivity of the coronal 
plasma in comparison to the photospheric one, the system of electric circuits 
is closed through the corona, if there is a magnetic connection between the 
two regions of the bipolar active region. Then a large scale DC electric field 
can establish in the corona which can be held responsible for the electron 
acceleration. 

Adopting plausible values of parameters in the solar atmosphere (as 
presented in this paper), an electric field of about — 2 mV /m is generated in 
the corona. Such a field can accelerate sufficient (depending on the plasma 
temperature more than 2 % of the total) electrons to kinetic energies of up 
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to ar ound 240 keV in typical time scales of less than a second (jAschwanden 



The resulting fluxes of energetic electrons are of the order of 10 36 U , 



as it is really observed with the RHESSI spacecraft (jWarmuth et al.l . 120071 ) . 

According to Ampere's law electric currents are the sources for magnetic 
fields. Basing on F e 10 36 V s and considering (as usually observed) two 
hard X-ray sources, the electrical current related with the X-ray emission 
can be estimated by I = eF e/2 = 8 x 10 16 A. Such a current localised in a 
coronal loop would be related to a magnetic flux density in the order of 
B « voi/is 1.44 kT, where no represents the vacuum permeability. On one 
hand these extremely high magnetic fluxes are not observed in the solar 
atmosphere, on the other hand the fluxes of energetic electrons deduced 
from the hard X-ray observations (e.g., RHESSI) are related with such high 
electric currents of the order of 10 17 A. In order to cancel the high magnetic 
fields induced in the corona by these high currents, the coronal currents 
need to be oppositely directed and cospatially located to each other. The 



system of electric circuits as discussed here (see Fig. 3(a)) has two electric 
currents (i.e., 15 and iq) through the corona. Note that the induced current 
is independent from L p h, i.e. the photospheric length scale (see Eq. ©). 
The important conclusions of the model presented in this paper are: 

1. The energy generated by the photospheric plasma motion is trans- 
ported electrically into the corona where it is transferred into the 
flare. 

2. The photospheric plasma motion has enough power to support a flare 
in the corona. 

3. The electric fields related with the potential drops at the coronal re- 
sistors are oppositely directed, so that the electrons are accelerated to 
two opposite directions leading to a double source of the hard X-ray 
radiation, as it is usually observed. 

The hard X-ray radiation is produced by br emsstrahlung of the energetic 
electrons in the dense chromospheric plasma ( Brown . 1971 . 19721 ). 



In summary, in the solar atmosphere electric circuits can be driven by an 
electric power supply induced by the photospheric motion in a bipolar active 
region. If these electric circuits are closed via the corona, an electric field is 
established along a coronal loop. Electrons can be accelerated along these 
large scale DC electric fields. The so produced energetic electrons can emit 
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hard X-ray radiation. Therefore the model presented is in good agreement 
with observations. 
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